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Abstract 



We perform infrared factorization of differential rates of radiative and semileptonic inclu- 
sive decays of the B meson in the end-point region of photon and charged lepton spectrum, 
respectively, in the leading heavy quark mass limit. We find that the differential rates are 
expressed in terms of hard, jet and soft functions, which satisfy evolution equations. Solving 
these equations, we find expressions for the moments of the differential rates in their end- 
point regions, which take into account all leading and nonleading logarithmic corrections in 
perturbation theory, as well as large nonperturbative power corrections. Expanded to the 
one-loop level, our predictions coincide with the results of existing lowest order calculations 
for B — > 'yXg and B — > WX U . Nonperturbative corrections appear in our formalism from 
the boundary value of the soft function in the evolution equation. The soft function is 
a universal process- independent function, which describes the distribution of the b quark 
in the B meson. Its behavior in the end-point region is governed by the nonperturbative 
asymptotics of a Wilson line expectation value. By considering the contributions of infrared 
renormalons, we find an ansatz for the Wilson line, which leads to a Gaussian model for 
the heavy quark distribution function. 
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1. Introduction 

Radiative and semileptonic decays of B mesons near maximum photon and charged lepton energy, 
respectively, are of special interest for the determination of electroweak parameters and the 
detection of new physics. They also present special theoretical challenges |] because of large 
perturbative and nonperturbative corrections from QCD. 

Recently, Bigi et a/.|J and Neubert || analyzed nonperturbative corrections to these processes 
in the A/M expansion, with M the heavy quark mass and A the usual QCD scale. They found 
that the leading term in this expansion is described by the same universal function, f(x), for both 
processes. This function is nonperturbative in origin, and describes the B meson decay outside 
the phase space available for the decay of a free b quark. They also related the first few moments 
of f(x) to hadronic matrix elements in heavy quark effective theory (HQET). As pointed out in 
||, 0, [5], however, perturbative corrections are also large in the end-point region x = E/E max ~ 1 
and need to be resumed. It was observed that the leading (double-logarithmic) corrections, 
(a s ln 2 (l — x)) k , exponentiate ||. Clearly, a complete picture of inclusive B decay in the end- 
point region requires a unified treatment of perturbative and nonperturbative contributions. Our 
goal in this paper is to develop such a formalism, based on well-known factorization theorems 
and evolution equations in perturbative QCD ||, and on the analysis of infrared renormalons 
0-0 in the evolution of Wilson lines, or path-ordered exponentials. 

In Section 2, we apply the factorization technique JTD], [II], |12[ to inclusive radiative decay 
B — > 7X5 in the photon end-point region. The relevant factorization is valid to leading power 
in 1/M, and involves two functions, one of which is the distribution f(x) describing soft gluon 
emission by the b quark, and the other a "jet" function J(x), describing collinear interactions 
of the outgoing s quark. The "universality" of f(x) is a direct consequence of factorization. In 
Sections 3 and 4 we show the close relation of the distribution and jet functions to Wilson lines. 
We then derive the evolution of these functions, which resums both leading and nonleading log- 



arithms as x —> 1. These equations are analogous to the usual GLAP [13] evolution equations, 
and follow directly from the renormalization properties of Wilson lines and quark fields. The so- 
lutions to these equations give an expression for dT-B-*-yx s /dx that includes all large perturbative 
corrections near the end-point. As in most evolution equations, they require a boundary value, 
which summarizes nonperturbative corrections. In Section 5, we show that the perturbative so- 
lution contains infrared renormalons, which introduce ambiguities in perturbation theory. These 
ambiguities occur precisely in the powers A 2 /[(l — x)M] 2 that we expect from nonperturbative 
effects. This analysis lead us to a nonperturbative estimate for the relevant Wilson lines, and as a 
consequence, to a Gaussian model for the nonperturbative behavior of the distribution function. 
Finally, we extend our analysis to the end-point region in semileptonic decays, and show that 
the only significant difference is that these decays require an extra integral over the phase space 
of the neutrino. 

2. Factorization in radiative decay of B-meson 

Let us consider the inclusive radiative decay B — > ^X s in the limit when the mass of the s quark 
is neglected. In the rest frame of the B meson we define a scaling variable x as the ratio of the 



photon energy to the mass of the b quark, 

The condition M\ > for the final state X s leads to the restriction [fj, |3|, 

< x < M B /M ~ 1.09, 

where Mb is the mass of the B meson. We are interested in the inclusive distribution, ^, 
with the energy of photon near the end-point, x ~ 1. The decay B — > ^X s occurs through 
the transition b — > sj, which is described by an effective hamiltonian [rj|. In the heavy quark 



limit, M — > oo, the momentum of the b quark is Q^ = Mv^ , where v^ is the velocity of the 
B meson. As was stressed in |2], [I], perturbation theory describes a free b quark decay in the 
region x < 1, and in order to penetrate inside the "window" 1 < x < Mb/M one has to take into 
account nonperturbative effects. In the rest frame of the B meson, the momenta of the photon 
and b quark, q^ and Q^, respectively, may be taken to have the following light-cone components 
(q± = 73 Oo ± 93) and q t = (q 1: q 2 )) 

M M 

q^ = (q + ,q-,q t ) = -j=(x,0,0) , Q, = y|(M,0). (1) 

Then, the momenta of the s quark, P = Q — q, is given by 

M , 

In the end-point region, x — »• 1, the s quark moves in the minus direction with a high energy 
~ M/2 and has a small invariant mass, P 2 = M 2 (l — x) <^C M 2 . Thus, for x ~ 1, the s quark can 
produce a jet of collinear particles accompanied by soft radiation, from which we expect large 
perturbative and nonperturbative corrections. This is exactly the same situation one encounters 
analyzing the behaviour of hadronic processes (DIS structure functions, Drell-Yan cross section) 
in perturbative QCD near the boundary of phase space HlQ| . That is why the resummation 



technique developed in Jul |TT], [L2j applies to the inclusive B decay in the end-point region. In 
particular, analyzing diagrams of fig.l contributing to the decay B — > '-/Xg we find the following 
three configurations of particle momenta associated with leading 1/M behavior (analog of leading 
twist in DIS): 

Hard(#) : k + ~ k_ ~ k t = 0{M) , 

Jet (J) : k + = 0(M(1 - x)) , k_ = O(M) , k t = O(M^T^) , 

Soft (5) : k + ~ A;_ ~ k t = 0(M(1 - x)) . 

For the x — > 1 limit, the hard subprocess gets contributions only from virtual particles. Notice 
that the minus and transverse components of the momenta of particles in the jet subprocess are 
much larger than those in the soft subprocess. That is, the jet subprocess carries almost all the 
P_ and Pt momentum, while the small P + momentum is distributed between the jet and soft 
subprocesses. In individual diagrams of fig.l, particles from jet and soft subprocess interact with 
each other. In the sum of all diagrams in the leading 1/M limit, however, the contribution of 
hard, jet and soft subprocesses may be factorized into the form || 10| 



1 dT M fP+ .,„„,„„ ,,„,_, „ M 



I \ dl + S(l+)J(P+ - l + )H(PJ) , P + = -Wl - x) . (2) 



1 j (IX V+ J{l+)min V 2 



Here, r 7 = T(b —> 37) = ^-^M 5 \VtbV^ s \ 2 Cj(M) is the partonic total width in the Born approx- 



imation |D| , and l + is the light-cone component of the total momentum of soft gluons emitted 
in the partonic subprocess. Let us make an important comment about the integration limits for 
1+ in (Q). The upper limit (l + ) m ax = P+ follows from the condition that the momentum of the 
s quark jet, P — I, be time-like. The lower limit (l+) m i n corresponds to the minimal energy of 
soft gluons in the final state. In perturbation theory soft gluons are emitted by b and s quarks 
into the final state and momentum conservation requires (l+) m i n = 0. In a realistic B meson, 
however, soft gluons may also take energy from the light components of the B meson. Thus, the 
minimal energy of soft gluons emitted by the b and s quarks in the partonic subprocess may even 
be negative nonperturbatively, 

(l + ) mi n = -(M B -M)/V2, (3) 

which leads a window for 1 < x < Mg/M. 

To separate the subprocesses H, J and S in momentum space one has to introduce || a 
factorization scale ft. The contribution of each subprocess depends on this scale, although the 
/i— dependence cancels in the differential rate @. In particular, the contribution of the hard 
subprocess, H, depends only on M and /i.Q 

3. The soft subprocess 

In emitting soft gluons, the b and s quarks behave as classical relativistic particles. That is, 
all effects of their interactions with soft gluons are factorized into eikonal phases given by the 



path-ordered exponentials [T(| [171 , or Wilson lines, P exp(i J c dz ■ A(z)), evaluated along their 
classical trajectories C with the gauge field A(z) describing the soft radiation. As a consequence, 
the contribution of the soft subprocess to the differential rate is given by a Fourier transformed 
expectation value of a Wilson line^ Q, [H| 

S(l+) = Mv + |~ ^t-jy-i+Wc((v ■ y)ii) , W c = (B\P exp (i ^ dz ■ A{z)^j \B) , (4) 

where the integration over y_ fixes the total momentum of soft gluons in the final state to be 
equal to l + and \B) denotes the B meson state. The integration path C in the definition of Wc is 
shown in fig. 2. It goes from — 00 to point along the quark velocity v, then along the light-cone 
minus direction to point y_, and then from y_ to —00 along —v. 

The following comments are in order. Notice that according to the definition, S(l + ) depends 
only on the properties of the bound state \B) and not on the particular short distance subprocess. 
This is why S(l + ) is a universal distribution. In fact, the function S(l + ) was previously defined in 



TT| in the analysis of large x behavior of the structure function of DIS. It can be easily checked 
that the heavy quark distribution function 0, [| coincides with the definition of S(l + ). Indeed, in 
the leading 1/M limit, heavy quark fields are equivalent to path-ordered exponentials ||19|| , and 
making this identification one gets S(l + ) = F(x) = f(k+) with l + = —Kx/\/2 = —k + /\^2 in the 
notation of refs. and , respectively. 

^n what follows the dependence of hard, soft and jet subprocesses on the coupling constant a s (n) is implied. 
2 Notice that the gauge fields are ordered along the integration path C and not according to time. However, 
on different parts of C path-ordering implies time or anti-time ordering pa] . 



In perturbation theory, the Wilson line Wc obeys the following relations fill , [18 

W c = W(y_v + p - ie) = W*{y_v + p - ie) = W{-y_v + p + ie) , 

which ensures the reality of S(l + ). Here, p is the factorization scale and the u —ie" prescription 
comes from the analogous property of the free gluon propagator and is thus of perturbative origin. 
Let us define a dimensionless variable z and rewrite (Q) as 

S(l + ) ee f(z, M/p) = r ^j^Wina/M -ie), l + = ^=(l-z). (5) 

J-oo Z7T \/2 



where a = y_v + M. Then, the u —ie" prescription immediately leads to the spectral property [11 

f{z,M/p)=0, forz>l, 

which implies that there is no "window", z > 1, in perturbation theory. For this window to 
appear one has to take into account nonperturbative corrections to the Wilson line Wc in (§). 

The heavy quark distribution function defined above, f(z,M ///), gets large perturbative 
corrections ~ ( ° g \-^ ) from the region 2: — > 1 corresponding to l + — > in (0). Analyzing the 
moments of f(z, Mj ' \x) in perturbation theory we find the following remarkable relation 



rM B /M 

f n (M/fi) = / dzz^fiz, M/fi) = W{-inn/M) + 0(l/n) , (6) 

Jo 

which means that the large n behavior of the moments of the heavy quark distribution function 
is given by the Wilson line expectation value Wc evaluated along the path C with the formal 
identification y-V+ = —in/M. This suggests that in the large n limit one has to treat Wc as a 
nonlocal functional of the gauge field [|Tl~1 , 18] , rather than to expand it into a divergent power 



series in y_ using the operator product expansion 0, [|. 

We now use the renormalization properties 0, |21], |22| of the Wilson line Wc as a function 
of (jl, which depend on the particular form of the path C. The integration path C of fig. 2 has 
two cusps at points and y and a light-like segment in between. As a consequence, the light-like 
Wilson line Wc obeys the renormalization group equation JTT], |l8l 

VWc = -{r cusp {a s )[log{p-ie)+log{-p + ie)} + r{a s )}W c , p = y_t; + /j e 7 ^ , (7) 

where T> = p-j^- + P(g)^-, 7 is the Euler constant and the anomalous dimensions T cusp (a s ) and 
r(a s ) are known to two-loop order |^1], [Tj|. Solving this equation, and using the relation (|), we 



find for the moments of the distribution function 



CpMno/n rjh \ h.n 

\ ^ 2T cusp (a s (k t ))\og^- + T(a s (k t )) 



Mm 



1 n 



where n = e 1 n and f^' has a nonperturbative origin, as the boundary value of the Wilson line 
in the solution of the RG equation (|7|). From this expression, we find that the moments obey the 



following evolution equation |TTj, [18 



V ,/„ ( \ f/p) = ( 2T cusp (a s ) log ^- + T(a 8 ) j /„ ( \ Up ) ( 9 ) 



where the r.h.s. depends on the normalization point through log 



Tlfl 

n M' 



4. The jet subprocess 

The jet subprocess J in eq.(||) describes the decay of the s quark with momentum P — I into a 
jet of collinear particles. The invariant mass of the jet is 2P_(P + — l + ) = M 2 (z — x) > 0, and 
the contribution of the jet subprocess depends only on this quantity and on the renormalization 
scale /i, 

J(P + - /+) = J(2P_(P + - /+), /i 2 ) = J(M 2 (z - x)^ 2 ) . (10) 

To take into account large perturbative corrections from the region z->iwe again consider the 
moments: 

JJM/u) =M 2 f dzz n - 1 J(M 2 (l - z),fi 2 ) . (11) 

JO 

After integration of the r.h.s., using the relation z n « e -n(i-z) - m ^q large n limit, the function 
J n depends only on two scales, M 2 /n and /i 2 , which allows us to write 

J n (M/fi) = V /M2 (l) + 0(l/n) . (12) 

Let us now substitute expressions fl5|) and (0) into (||), using @, 
1 dT r M B /M 

— — = / dz f(z,M/fi) M 2 J{M 2 {z-x)^ 2 ) H(M/fi), (13) 

1 <y CLX JX 

and consider the moments of the differential rate defined as 

1 i-Mb/M ^F 

M n (B -h. 1 X S ) = — dx x^ 1 — . (14) 

T 7 Jo dx 

The behavior of dT /dx in the end-point region x ~ 1 corresponds to the large n asymptotics of 
the moments -M„(B — > jX s ). We notice that for z, x ~ 1 one can replace J(M 2 (z — x),/i 2 ) ~ 
J(M 2 (1 — x/z),fi 2 ) in (|T3|) . Then, the moments of the differential rate factorize [[L(J into the 
product of moments of the distribution function /„ and the moments of the jet J n defined in (Bf) 



and (11 



M n (B -> 7 X S ) = f n {M/n)J n {M/n)H{M/n) + C(l/n) (15) 

The condition that the l.h.s. of this relation does not depend on the renormalization point /x 
leads to the following RG equations, 

VJ n {M/^i) = \2T cusp {a s ) log ^--2 7 (a s ) J J n (M/fj,), (16) 

VH(M/y) = {-2T cusp {a s )\og^ + 2 1 {a s ) + T{a s )^H{M/^, 



where the functional form of J n , eq . (|12|) , and independence of H on n were used. The anomalous 



dimensions entering into these equations have the following one-loop expressions [II, |18 

F C usp(a s ) = —C F , T(a s ) = S -C F , j(a s ) = --—C F . (17) 

TV TV 4 7T 

Here, 7(« s ) coincides with the one-loop quark anomalous dimension in the axial gauge. To 
understand this property we notice that the collinear subprocess J can be defined as a cut 
propagator of the s quark in the light-like axial gauge (n ■ A) = with n M = (0, y_, 0). Indeed, in 



this gauge the interaction between soft gluons and particles from the jet subprocess is supressed 
in each diagram of fig.l and the factorization (|2|) is manifest. Then, the evolution equation ( |16[) 
describes the renormalization properties of a cut quark propagator in the light-like axial gauge. 
Solving the RG equations (H) and flilf ) for S n , J n and H, we find after some algebra our 
complete expression for the moments of the decay distribution, 




/f C,(a s (M)) 

2 / ^ W« s (A; t )) + T(a s (My)) + i(a 8 (My/y)) 

JMv kt. 



dy 

n /n 



I My k t 

where C 7 = 1 + 0(a s (M)) takes into account log (n) corrections to M. 

5. Infrared renormalons in Wilson lines 



For large n, the integrals in ([18]) get dominant contribution from the region y ~ 1/n, and for 
n ~ M/A we encounter the singularities of the coupling constant. This means that perturbation 
theory becomes ill-defined in the end-point region of the photon spectrum, and we have to take 
into account nonperturbative corrections to dTB-,-yx s /dx. As discussed in J7|, both perturbative 
and nonperturbative contributions to physical quantities are ambiguous at the level of power 
corrections and only their sum is unique JF], |23|. This allows us to understand the structure of 
nonperturbative effects by analyzing the ambiguity of the perturbation series associated with the 
so-called infrared renormalons. 

In the end-point region we expect large nonperturbative corrections to the Wilson line entering 
the definition (^) of the distribution function. Let us examine the contribution of IR renormalons 
to Wc by performing an "improved" calculation of Wc in perturbation theory. Using the defi- 
nition (|j) we find the one-loop contribution to Wc [IS] and use the nonabelian exponentiation 
theorem for path ordered exponentials ||24] to get 



Wo = exp CV-/ " ^MUt- £A (1 - e«)(l - e"«) , (19) 



where y = (0+, y_, 0), d D k = dk + dk_d D ~ 2 k, and where \i is the scale parameter of the dimensional 

regularization with D = 4 — 2s. The choice of the argument of the coupling constant a s is 

determined by higher order corrections to Wc to be the transverse momentum of the gluons 

17J |25j . Now let us substitute a s = a s (k 2 ) = / °° da exp (— a (3q log j?) into fli9|) and perform the 



integration over gluon momenta to get 

(4tt^ 2 /A 2 ^ 



exp CV l F 7/ { f^ da (Ay_v + ) 2 ^°T(-2a{3 ) cot(7ra/Sb)(l - a(3 ) . (20) 

\ r(l-e) Je/0o 



W c 

/A 



Because (Ay_v + ) 2crl3 ° = exp(—a/a s (l/y^.v + )), the exponent has a form of an inverse Borel trans- 
formation 0-0. Integration over small Borel parameter a corresponds to large transverse mo- 
menta k 2 and gives rise to ultraviolet poles in s. After their subtraction in the MS scheme, one 
verifies that Wc satisfies the RG equation (|7|). We now recognize that away from a = s/j3o the 
integrand in fl20"|) has infrared renormalon singularities generated by the T— function at points 



a* = 27f, 4-, 2§~, •••• Thus, to give a meaning to the perturbative expansion one has to fix a 



prescription for integrating these singularities |J. Different choices of the prescription lead to 
results which differ in power corrections of the form (Ay_) 2cr *^°. The leading power corrections 
correspond to the left-most singularity. Note that the singularity of the T— function at a* = ^4- 
is compensated by the cotangent, so that the leading IR renormalon appears at a* = «-. Thus, 
resummed perturbation theory generates (but fails to describe uniquely) 0(A 2 y 2 _ ) power correc- 
tions to the Wilson line Wc- At the same time, this means that nonperturbative effects should 
also contribute to 0(A 2 y 2 _ ) power corrections to make the Wilson line 

W c = W pert . (y_ ) W nonpert . (j/_ ) 

well defined. In summary, the exponentiation of perturbative corrections to W per t.(y-) implies 
that IR renormalons appear in the exponent of W per t.(y-)- Hence, in order to compensate their 
contribution one can choose the "minimal" ansatz for nonperturbative part W nonper t.(y-)'- 

W nonpert {y-) = exp (-(y_v + f/l 2 + 0{yl)) , (21) 

where I is a nonperturbative correlation length. Notice that the expansion in the exponent does 
not contain a term linear in y_. 

Let us test the ansatz of eq. (pl|) by performing an expansion of the Wilson line Wc defined 
by (U) in powers of y_. We denote a straight Wilson line as <& v (x) = P exp(i f®^ dsv ■ A(vs + x)) 
and apply the identity Pexp(z/ dsy_A + (ys)) = exp(—y_D + ) exp(y_<9 + ) to get 

Wc = (B\& v {y_)Pe l h ^-^(^$,(0)15) = 1 + E -y-(B\^M {*D + ) k $,(0)|£?> 

fe=i 

where D^ = <9 M — iA^ is the covariant derivative in the quark representation. According to its 
definition, the Wilson line satisfies the relation (v ■ D)<& v (x) = 0, which is similar to the equation 
of motion of the heavy quarks in HQET |26] . This allows us to apply HQET machinery for the 
evaluation of coefficients in the expansion of Wc- 

W c = l + oi(*y_u+) + a 2 {iy-v + ) 2 + a 3 {iy-V+) 3 + ... , (22) 

where 

oi = v- + l {B\& v {iD + )§ v \B) = , a 2 = -\(B\& v {iD) 2 $> v \B) , a 3 = —(Bl^F^D^JB). 

o lo 

We recognize that our ansatz (F2TJ) is consistent at 0(y 2 _) with the small y_ expansion of Wc- 
Note that in the leading 1/M limit, the heavy quark fields h v are proportional to the Wilson line 
<& v , which is why the coefficients a n are equal to analogous fundamental parameters in HQET 



0, Hj. In particular, we can identify the correlation length I entering into ( 21 ) as 

I 2 = -o 2 = -6/Ai = 6//4 • (23) 

where the value of /i 2 = 0.52 ± 0.12(GeV) 2 was estimated using QCD sum rules in p7| (which 
implies a particular choice for the perturbation series). Substituting the nonperturbative estimate 
(21) for the Wilson line Wc into (J5|), we find a distribution function f(°'(x), which can be used 



as a nonperturbative input for the evolution equations (g). If we use only the small y_ expansion 
(|22| ) we get an expression for p°'(x) as a series of derivatives of 5(1 — x)— function, equivalent 

8 



to that proposed in 0, Q , which requires us to guess the shape of the function in the end-point 
region x ~ 1. On the other hand, integrating the nonperturbative ansatz ([H]), inspired by the 
IR renormalon analysis, we find the function Q 

which describes a Gaussian distribution around x — 1 with the width a. Notice however, that 
since we neglected 0(y^_) corrections to the exponent of Wc in Q2T1) , this function does not vanish 



outside the physical region x > Mb/M and (p4|) is valid only at the vicinity of x = 1. To get 



a better approximation for the distribution function one needs a more realistic estimate for the 
Wilson line expectation value. Once we know the distribution function f(°>(x), we evaluate its 
moments, f^ = f dx x n ~ 1 f^°\x), and substitute them into (||) and (|18|). 

6. Factorization in semileptonic decay of the B meson 

The factorization for B — > jX s was based on the kinematical requirement that, near the end 
point of the photon spectrum, the s quark must form a jet with large energy and small invariant 
mass. The same situation appears in semileptonic decay B — > lvX u in the end-point region 
of the charged lepton spectrum. In this case, the decay occurs through an effective local vertex 
ujfj,(l — 75)6- ^7 M (1 — 75)^- Let I and v be the momenta of the leptons, and Q and P the momenta 
of the b and u quarks, respectively. Standard notations for the scaling variables in the rest frame 
of the B meson are, 

2E l W 2 2W 

where W = l + v = Q — q is the momentum transferred to the leptons. Phase space for these 
variables is 

< x < x m = M B /M, < y < xx m , y/x + x < y < y/x m + x m . (26) 

In analogy to the radiative decay, we examine the energy and invariant mass of the jet created 
by the u quark, 

P 2 = M 2 (l-y + y), P = ^-(2-y ). (27) 



From the definition (|26f) of phase space we find that 

P 2 < M 2 (l - x)(l - y/x), p <M-(l- y ) + M-(l - x)(l - y/x) . (28) 

We conclude that near the end point of the lepton spectrum, x — > 1, with y < 1, we meet the 
same kinematical situation as in the B — > ^X s decay. Thus, once again, in the leading 1/M limit 
the dominant contribution to the differential rate of the decay comes from diagrams containing 
hard, jet and soft subprocesses. As in the previous case, the contributions of these subprocesses 



3 The same function has been proposed in 



are factorized from the partonic subprocess o"o ~ (Q ■ v)(l ■ P) = ^r~(x — y)(yo — x), and for the 
triple differential rate we get an expression similar to (fj), 

1 d 3 T M r p + 

-— — = (x- y)(y Q - x) - I dl + S(l + )J(P + - l + )H(P.) . (29) 

Ti dxdydyo v + J(i + ) min 



G 



where Ti = j^\V u b\ 2 M 5 , and (l + ) min was defined in (Q). Note that the prefactor in ( |2"9"D sup- 
presses the "dangerous" region y — ► x ~ 1. In this region, P is forced to zero by ( p8|) and 
factorization fails because the outgoing quark has vanishing energy. We note that relations (0) 
and ( p9[) were found in a frame, where the outgoing quark has momentum P + <C P_. For the 
radiative decay, this frame is fixed by the particular choice (H) of the photon momentum q. For 
a semileptonic decay, we use the same choice for the lepton momentum, I = -^(x, 0,0). Let us 
express P + and P in terms of the scaling variables (E5|). Using the relation P + <C P_ we make 
the following approximation: P_ ps P_ + P + = \^2P = y/2(M — Wo) and get 

M P 2 M 1-v 

P. = ^(2-«,), P + = — = -^(1-0, { = 5-^, (30) 

where the scaling variable £ is an analog of x in the radiative decay, eq.(§). However, in contrast 
to the previous case, £ depends on y and y 0) an d to calculate the differential rate dT B ^i PXu /dx 
one has to integrate (p9|) with respect to £ in the physical region x < £ < Mb/M. Using the 
kinematic relations (|30| ) and the definitions ([|) and flip]) of the subprocesses we can rewrite (p9"D 
as 

i *r = (I _ , )(yo _ x) r-<- dz f L K ±y. HMH2 _ m)(z _ «,, v (m^m) . 

T; dxdydyo h \ (i> J \ H J 

(31) 
Notice that if we replace subprocesses by their lowest order expressions H = 1, / = p°^(z), 
J = <5(M 2 (2 — yo)(z — £)), we get a result which contains all nonperturbative corrections, and 
which coincides with analogous expression in 0. 

In order to analyze large perturbative corrections to dT /dx in the end-point region, we change 
variables from y to £, and integrate Q3"T| ) with respect to £ and y , neglecting terms that vanish 
as x *■ i , 

1 rlT f 2 rM B /M i-Mb/M 

Tdx' = I dy ° (2 " yo)2{m ~ 1} I di k ^ (32) 

x/(^)m^(2-, )( 2 -(),^)i(^) + C)(1-,). 

Comparing this expression with ( |13| ) we notice that — ^(t^tF) is similar, apart from y — integral, 



to the differential rate f-% of the radiative decay, eg. (|T3|) . This suggests the appropriate mo- 
ments for semileptonic decay in analogy with (0), 

i-Mb/M d / ^ dV\ n fM b IM AT 

M,(B - IP*,) = - / <^U ^T"- (33) 

jo dx \l i dx J L i Jo dx 

As in the case of radiative decays, we evaluate the moments of 4^, eg. (|32|) , and find in the large 
n limit the following relation 

M(1 ; g - ) ) +0 (l/n)(34) 



M n (B - /z/X u ) = /J— ) jT'dar, ^(1 - *„) jJ MV1 ^ )p[ 
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where we have changed variables to y$ = l + x u « 1 + ^ff with E u the energy of outgoing neutrino 
in the rest frame of the B meson. Here, the moments f n and J n are identical to whose in the 
radiative decay and satisfy the evolution equations (|D and fllB]). The only difference with ( |T5D 



is that one has to integrate in (p4|) with respect to the energy of the neutrino, x u , and take into 
account the x u — dependence of the collinear and hard subprocesses. This additional integration 
does not affect leading double logarithmic corrections (a s \og 2 n) k , but it does affect nonleading 
terms. Note that in fl34|) the prefactor suppresses the end-points x v = and x v = 1, which 
correspond to the limits of soft neutrino and outgoing quark, respectively. Thus, our factorized 
expression ( |3~3| ) takes care by itself near the "dangerous" point x v = 1, where IR factorization 
fails. 

Let us compare our predictions flip] ) and fl3~3| ) for resummed large perturbative corrections with 



the results of one-loop calculations |28|, [15| of the differential rate for B — > ^yX s and B — > lvX u in 
the end-point region. Solving the evolution equations (|]) and (0) we find one-loop expressions 
for the subprocesses: 

MM/,) = i + S^f-ltfJEL + log "" 




Mn Mn c 

W') = 1 + ^|2.og^ + |.og-^) (35) 



H[Km = 1 + ^(-log^-flog-£) 



where we have omitted constant terms and ffl. Substituting these relations into ([TED and (H), 
we obtain the one-loop expression for the moments of the differential rates 

M n ~ 1 + — C F (- \og 2 n + A \ogn + const.) (36) 

where coefficients in front of nonleading log n— term are 

Ab^ 7 x 3 = - , A B ^ WXu = — (37) 



in accordance with the one-loop results of ref . [Q and [|15j . 

Using the evolution equations @ and (|16D , one can represent the expression ( p^) in a form sim- 
ilar to (|18|). It is more interesting however, to consider the ratio of the moments .M n (B — > ~fX s ) 
and A4 n (B — > WX U ) defined in (^) and (|34l) , respectively. We find that the moments of the 
heavy quark distribution function, /„, cancel in the ratio and we get 

M n (B -> WX U ) Q /-i ( fi dy [My dk \ 

— -— tty = TT \ dx u x u {l - x u ) exp 2/ — / —T cusp (a s (k t ))\ 38 




7 (a s (M^/yn /n) ) + 2 7 (a s (My)) + r(a,(Mj/)) 



x exp 

where Cj = 1 + 0(a a (M)) is an analog of C 7 , defined in (0), for the case of semileptonic decay. 
Thus, all nonperturbative corrections cancel in the ratio of the moments of the differential rates 
and this allows us to calculate (p8|) perturbatively. We may then compare it with the ratio of 
experimental data using the definitions of the moments, (|H|) and (|33|), and isolating the ratios of 
the elements of the Cabibbo-Kobayashi-Maskawa matrix contained in the prefactors T 7 and IV 
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7. Conclusions 

In this paper we performed infrared factorization on the differential rates of radiative and semilep- 
tonic inclusive decays of the B meson in the end-point regions of the photon and the charged 
lepton spectra. We found that in the leading 1/M limit the differential rates are expressed in 
terms of hard (H), jet (J) and soft (S) functions which satisfy evolution equations. 

Solving the evolution equations we found expressions for the moments of the differential 
rates in the end-point region, fllBp and (|3"8"|), which take into account all leading and nonlead- 
ing logarithmic logn corrections in perturbation theory, as well as large nonperturbative power 
corrections in the leading 1/M limit. Expanding these expressions in powers of the coupling 
constant and using one-loop results (|17]) for the anomalous dimensions entering the evolution 
equations, we have shown that our predictions coincide to the lowest order with the results of 
one-loop calculations for both processes. 

Nonperturbative corrections appear in our formalism from the boundary value of the soft func- 
tion /(°)(x) in the evolution equation (Q. The soft function is the universal process-independent 
function which describes the distribution of the b quark in the B meson. We established that the 
behavior of f^(x) in the end-point region x ~ 1 is governed by the nonperturbative asymptotics 
of the Wilson line expectation values. Considering the contribution of infrared renormalons, we 
found a nonperturbative ansatz for the Wilson line which led to a Gaussian model fl24 ) for the 
heavy quark distribution function f(°\ 
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Fig.l: Unitarity diagram contributing to the differential rate of the radiative 
inclusive B — > ^X s decay in the leading 1/M limit in the axial gauge. We 
use solid lines for quarks, dotted lines for gluons, wave lines for photons and a 
dashed line for the final state. 



j/ = (o+,j/_,5) 
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Fig. 2: Integration path C entering into the definition (f|) of the Wilson line 

W c . 
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